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In this paper flow patterns of the labyrinth seal are experimentally 
investigated for making a mathematical model of labyrinth seal and to obtain the flow 
induced force of the seal* 

First, the flow patterns in the labyrinth chamber are studied on the 
circumferential flow using bubble and on the cross section of the seal chamber using 
aluminum powder as tracers. And next, the fluid force and its phase angle are 
obtained from the measured pressure distribution in the chamber and the fluid force 
coefficients are derived from the fluid force and the phase angle. Those are similar 
to the expression of oil film coefficients. 

As the results it is found that the vortexes exist in the labyrinth chambers and 
its center moves up and down periodically. The pressure drop is biggest in the first 
stage of chambers and next in the last stage of chambers. 


INTRODUCTION 

Labyrinth seal are used as the non-contact type seals between rotor and stator 
of fluid machineries, such as, steam turbines, gas turbines, compressors and so on. 
To improve the efficiency of these fluid machineries, pressure of the working fluid 
and rotating speed of rotor tend to become high while the clearance of labyrinth seal 
is designed narrow to prevent for the leakage flow. As a risk of seeking such a high 
performance machinery, unstable vibrations of rotor occur and they draw an industrial 
important problem. The labyrinth seal is considered as one of these vibrations. 
Unstable vibration induced by the labyrinth seal was first pointed out by 
Alford(ref. 1), and after that it has been studied by many investigators. Among them, 
Boyman & Suter(ref . 2) made experimental studies and sketched the flow patterns. In 
this study, they had interested in the static characteristics of the labyrinth seal, 
but not in the dynamic characteristic, so, they did not concern with the flow induced 
force of labyrinth seal. On the other hand, Benckert & Wachter(ref. 3) had interested 
in the dynamic characteristic and they measured the hydraulic spring coefficients by 
using the static method. That is, rotor was fixed eccentrically toward casing and the 
hydraulic spring coefficient was measured by static way. The weak point of this study 
is not to be able to consider the effect of damping and its effect can not be 
neglected. On the contrary, Wright (ref. 4) investigated the instability condition of 
rotors induced by the labyrinth seal force for two seal strips model under many test 
conditions and also measured the destabilizing fluid force which included the damping 
coefficient. This study is valuable because it considers the hydraulic spring and 
damping coefficients, however it lacks generality because the rotor has only two seal 
strips . 
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The purpose of this paper is to study the flow patterns, which is needed 
mathematically to model the flow in the labyrinth seal. Namely, the flow pattern of 
the cross section and the circumference of the labyrinth chamber are observed. Bubble 
is used as a tracer for the visualization of the circumference and aluminum powder is 
also used in that of the cross section, and the photographs of them are taken. 
Moreover, the flow induced force is measured under the condition of both spinning and 
whirling motions of the rotor and measuring the change of pressure and of the phase 
angle in the labyrinth chamber. The eight coefficients, which represent the flow 
induced force of labyrinth seal by the function of velocity and displacement, i.e. 
they are correspond to the coefficients of the journal bearing, are obtained. 


TEST APPARATUS 

The test apparatus to observe the flow patterns and to measure the fluid force 
is shown in figure 1 . In this figure the rotor is driven by the three-phase motor 
@ through the variable-speed drive for spinning motion of (2) • The rotor is 
supported at the upper and the lower part by the eccentric bearings (4) • The 
dimensions and directions of both eccentricities are set to be equal. Therefore, the 
rotor is whirled by rotation of the eccentric bearing. The eccentric bearings are 
driven by the three-phase motor (Q through the variable-speed drive (Jl) for 
whirling motion. Since the rotor and the eccentric bearing are driven by different 
motors, the rotor can be spun in the same direction as whirl direction (forward) and 
in the opposite direction (backward), and these spin and whirl speeds are 
continuously changed from 1,6(Hz) to 5. 2 (Hz) . 

The schematic piping of the hydraulic circulating system is shown in figure 2. 
The working fluid is compressed and transported from the tank 1 to the tank 2. The 
tank 2 acts as an accumulator, and is installed to lower the pulsation of pressure by 
the pump using the air which is in the upper part of the tank 2. The working fluid is 
transported from the tank 2 into the test casing through four pipes, and is 
transported to the measuring part through the rectifier. The test seals have two or 
six strips and are the straight-through type seals. After that, the working fluid is 
returned to the tank 1 by natural draining. Two test casings are prepared, i.e. one 
is made of acrylic resin for the observation of streamline and the other is made of 
stainless steel for the measurement of the flow induced force. 

The rotating speed of the rotor is measured by the pulse of the notched disk 
which is installed on the upper end of the rotor. The whirl frequency of the rotor is 
measured by the measurement of the whirl amplitude of the rotor. The change of 
pressure by the whirl of rotor is measured through the needle of injector, which is 
fixed in the measuring point of the test casing where the middle point between two 
confronted strips corresponds, by the semiconductive pressure transducer. The 
measured signal current is sent to the signal analyser through the D.C. amplifire and 
is analyzed. 


EXPERIMENT I (Observation of flow patterns) 

Two kinds of experiment, that is, experiment 1-1 and experiment 1-2 , are done 
using the acrylic resin test casing in experiment I. And two kinds of test seal are 
prepared, i.e. one has two seal strips and the other has six seal strips. The flow 
pattern of the cross section of the labyrinth chamber is observed in experiment 1-1 
and the circumferential flow pattern is observed in experiment 1-2 . From the limit 
of stress of the test casing, two kinds of entrance pressure, i.e. 122.5(kPa) and 
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147.0(kPa), and three kinds of spinning speed O(Hz), 2(Hz), and 3(Hz), are tested 


Test procedure of experiment 1-1 (Cross section) 

To observe the cross section of the labyrinth chamber, the stroboscopic light is 
thrown through the slit whose clearance is 1.5 (mm) as shown in figure 3, and the 
reflected light of the tracer of aluminum powder is photographed. 


Test procedure of experiment 1-2 (Circumference) 

To observe the circumferential flow of the labyrinth chamber, the light of 
reflector lamp is used and the motion of vortex in the labyrinth chamber is 
photographed using bubble tracer. 


EXPERIMENTAL RESULTS AND CONSIDERATION 

Figures 4 and 5 show the flow patterns of one chamber seal, in which the flow in 
cross section is shown in figure 4 and the circumferential flow is shown in figure 5* 
The flow which through the first seal forms the vortex row in the labyrinth chamber. 
The center of vortex row does not stay and moves sinusoidally as shown in figure 5> 
and the sinusoidal wave moves in the same direction as spin, direction at a lower 
speed than rotating speed of rotor as shown in figure 6(a). This phenomenon is clear 
from three photographs of figure 4 that the cross sections of labyrinth chamber are 
taken. It is found from these photographs that the vortex moves up and down against 
flow as shown in figure 6(b). 

Figures 7 and 8 show the flow pattern of 5 chambers seal. The flow in cross 
section is shown in figure 7 and the circumferential flow is shown in figure 8. As 
the case of one chamber seal, the flow in each chamber also forms sinusoidal vortex 
row. The periods of vortex rows in 5 chambers are almost the same, but there is no 
particular regularity between the phase angles of the vortex rows. Bubble which 
through the first seal is big. It means that the pressure drop in the first chamber 
is high. This result corresponds to the measurement results of the pressure drop in 
each chamber. As the case of one chamber seal, it is found from figure 7 of the cross 
sectional photograph that there are vortexes. 


EXPERIMENT II (Measurement of fluid force) 

Two kinds of experiment, that is, experiment II-1 and experiment II-2, are 
performed using one chamber seal and 5 chambers seal in experiment II. The fluid 
force F which acts on the rotor, and the phase angle <J) between whirl and fluid force, 
are measured in experiment II-1 . The spring coefficients and the damping 
coefficients of the working fluid are measured in experiment II-2. 


Experiment II-1 (Fluid force and phase angle) 

The rotating speed is set to the whirl speed, and the change of pressure and its 
phase angle to the whirling motion are measured in each chamber. In case of one 
chamber seal, three points of figure 9(b), that is,(T),(2), and (5), are used and in 
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case of 5 chambers seal, seven points of figure 9(c), that is, CD to (7) are used. 
And the position of these measuring points is (III) of figure 9(a). The changes of 
pressure in each chamber are measured from ( 5 ) to (o) of figure 9(c) in order. In the 
analysis by the signal analyser, the 10 Hz-low pass filter is used to protect the 
noise of high frequency. Parameters are set to be entrance pressure, eccentricity, 
rotating speed, and direction of rotation. Conditions of experiment II-1 are shown 
in table 1 . 


Experiment II-2 (Spring and damping coefficients) 

The rotor is spun at the fixed eccentricity without whirling motion, and the 
circumferential pressure distribution is measured. As in the case of experiment II-1 , 
the measuring points are (T) to © in the position of (III) in figure 9» So, the 
directions of static eccentricity are changed one after another, its procedure is 
(III)-*-(IV)-»-(I)-KlI) in figure 9(a) > and the pressure of each eccentric state is 
measured at the measuring points. Consequently, the circumferential pressure 
distribution of the statically eccentric rotor is measured. The procedure of 
measurement is that, first, the pressure in each chamber is measured one by one when 
the rotor's eccentric direction is one of four directions ((I) to (IV) in figure 
9(a)) and next the direction of eccentricity is changed and the pressure is measured. 
Conditions of the experiment II-2 are shown In table 2. 


ANALYSIS OF MEASURED DATA 

The data obtained in experiment II are the maximum change of pressure AP max , the 
phase angle of pressure to the whirl motion /_§ and the change of the circumferential 
pressure distribution AP' when the rotor is spun without whirl. So, first, the x and 
the y directional force F xo and F yo are obtained from AP max and [<$> when the rotor is 
whirling. Next, the principal spring coefficient and the cross coupled spring 
coefficient Ky X are obtained from AP'. And last, the principal damping coefficient 
Cyy and the cross coupled damping coefficient C X y are calculated using the equation 
of fluid force and the expression is similar to that of the journal bearing. 

It is assumed that the circumferential pressure distribution in the whirling 
motion is circular and the center of it is Op as shown in figure 10. Thus, the change 
of the circumferential pressure distribution in the chamber AP is given as; 

AP = AP mcix C0S(y ♦ P -Ji) (1) 

The x and y direction components F xo and Fy 0 of the fluid force are obtained by using 
equation (1), 

tilt 

F xo =-/ q (AP COS ¥ )iR s d f 

= -iR s 7rAP max COS(7C-f) (2) 

F yo = -/ 2 o 7t (APSIN f )j?R s d? 

= -JtR s 7LAP max SIN(7C-tf> ) (3) 
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Then the fluid force and its phase angle acting on the rotor are as follows; 



where R s is radius of rotor and X is strip pitch. 

On the other hand, and Ky X are obtained from the change of the 

circumferential pressure distribution Ar f when the rotor is spun and no whirl with 
the fixed eccentricity, 

K X x = i-Jo 7l ( AP 'COS?)JJR s d ( f (6) 

1 ( 2K 

K yx =^rJ o (AP'SiN?)JtR s d<f (7) 

where e is eccentricity of rotor* AP 1 are measured for four cases, that is, ^=0, f- 
tt/ 2, y=7r, and f=3'n'/2, and are approximated to a circular distribution, as the case of 
AP, as follows; 

AP'=AP' max COS(<f + <!>’- K) ( 8 ) 

where AP f max is the maximum change of pressure and (J> ! is the phase angle between the 
position of AP * max and the maximum clearance. 

So, the principal spring coefficient K xx and the cross coupled spring coefficient Ky X 


are obtained as follows; 

Kxx'i^s^P'maxCOStTT-*) (9) 

Kyx = -^JlR s 7CAP , mQ xSIN(n:- £) do) 

Last, the expressions of fluid force induced by labyrinth seal are obtained as 
follows; 

F x = “(K xx * ^n^xy ~ ^n^xx 

V - < Kyx* " n c yy > * - <K yy- "nSx > y (12) 

where F x and Fy are the function of time and (% is critical speed of rotor system. 
Then equations (11) and (12) are considered at the same time as equations (2) to (5) 
as follows; 

F xo- ~ K xx e “ W n^xy e ^ 3 ) 

F yo = “^yx e ” w n^yy e ^4) 


Substituting equations (2), (3), (9) and (10) into equations (13) and (14) , the cross 
coupled damping coefficient c xy and the principal damping coefficient Cyy are 
obtained. 
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EXPERIMENTAL RESULTS AND DISCUSSION 


Experiment II-1 (Fluid force and phase angle) 

Figure 11 shows an example of wave mode of the signal analyser* The upper part 
is the wave of the rotor's whirl, and the lower part is that of the pressure change 
in the chamber. The noise of this wave is pressure disturbance induced by turbulent 
flow. It is found from figure 11 that there is phase lag between whirl and fluid 
force. 

Experimental results of one chamber seal are shown in figures 12 and 13* In the 
case of one chamber seal, the fluid force increases with the increase of eccentricity 
in every condition of both entrance pressures and whirl directions. In many cases, 
the fluid force increases with the increase of rotating speed. But in the cases of 
the high entrance pressure and forward precession the fluid force tends to decrease 
with more than a certain rotating speed. The phase angle <J) is constant, of which 
value is about 3Tr/2(rad) in all cases. 

Experimental results of 5 chambers seal are shown in figures 14 to 19* First, 
the discussions about the forward precession are as follows. The fluid force in each 
chamber increases with the increase of the entrance pressure. The fluid force in the 
1st chamber is more remarkably affected than the others by the rotating speed. The 
fluid force in the other chambers varies hardly or increases slightly with the 
increase of the rotating speed. And the fluid force increases with the increase of 
the eccentricity. The change of fluid force is shown in figure 16. This is the case 
that the rotating speed is 3*1 (Hz) and the tendencies of the fluid force do not 
change in the other cases. The fluid force in the 2nd chamber decreases, and in the 
other chambers after the 3rd chamber it gradually increases. The fluid force of the 
whole labyrinth seal F increases with the increase of rotating speed, entrance 
pressure or eccentricity (fig. 15)* The phase angle <f> is constant at about 3?r/2(rad) 
and has no relations to entrance pressure, rotating speed, and eccentricity, except 
for the case in the 2nd chamber. The phase angle in the 2nd chamber is changed by 
conditions. The reason of this phenomenon is considered as follows. Since the fluid 
force in the 2nd chamber is small, there is no remarkable difference between the 
pressure disturbance induced by turbulent flow and the change of pressure induced by 
whirl. Then, to obtain the accurate phase angle is difficult. In the case that 
eccentricity is large i.e. 1 .0(mm) and the fluid force is greater than the other 
cases, the phase angle ([> is constant by about 3'n'/2(rad). Thus, it is considered that 
the phase angles 4) in all chambers are constant by about 3 , n'/2(rad) and have no 
relation to parameters. In case of backward precession, the tendency of the fluid 
force is almost same as the case of forward precession. And the fluid force of 
backward precession is more remarkably affected than that of forward precession by 
the rotating speed. 


Experiment II-2 (Spring and damping coefficients) 

Part of experimental results of the spring and damping coefficients are shown in 
figures 20 to 31 * Variation in each figure may depend on the accuracy of instrument 
and measurement. and Ky X have no relation to the rotating speed of rotor. 

increases with the increase of entrance pressure, but Ky X has no relation to entrance 
pressure. K-^- and Ky X decrease during the working fluid flows from the 1st chamber to 
the 5th chamber. Cyy and C X y have no relation to entrance pressure and rotating 
speed. And C X y increases during the working fluid flows from the 1st chamber to the 
5th chamber. 
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CONCLUSIONS 


In this paper, the flow patterns of one chamber and 5 chambers seals are 

investigated, and the fluid force and its phase angle, and the coefficients of fluid 

force are experimentally obtained. The conclusions are as follows: 

(1) In the chamber, the vortex forms the sinusoidal wave and its wave moves in the 
same direction of rotation in the circumferential direction. Observing a 
certain cross section of seal chamber, the center of vortex moves up and 
down periodically. 

(2) The fluid force and its phase angle of one chamber seal are qualitatively agree 
with 5 chambers seal. 

(3) The fluid force increases with the increase of entrance pressure, rotating 
speed or eccentricity. 

(4) The fluid force in the 1st chamber and the 5th chamber of the 5 chambers seal, 
that is, the fluid force of entrance and exit part of labyrinth seal, is great, 
and the fluid force in the middle region of labyrinth seal is small. 

(5) The phase angle of fluid force is constant and has no relation to all 
parameters . 

(6) The principal spring coefficient increases with the increase of entrance 
pressure, and it decreases from the 1st chamber to the 5th chamber. The damping 
coefficients have no relation between entrance pressure and rotating speed. 

This study is supported by Scientific Research Fund of the Ministry of Education 

No. 56550180 (1981-1982) and No. 585501 72 (1983). 
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TABLE 1. - EXPERIMENT II -1 
EXPERIMENTAL CONDITION 


TABLE 2. - EXPERIMENT ft -2 
EXPERIMENTAL CONDITION 


| ENTRANCE PRESSURE tkPa ) 








122. S 147.0 171. S 196.0 












ENTRANCE PRESSURE (kPa) 


122.S 147.0 171.5 196.0 
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Figure 2. - The schematic piping of the circulating system 
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Figure 3. - The schematic layout of' the visualizational apparatus. 
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Figure 4. - Cross section of one chamber seal. 



Figure 5. - The movement of vortex row in one chamber seal- 
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of 5 chambers seal. 
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Figure 8. - The movement of vortex row in 5 chambers seal 



AMPLITUDE 
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O' : CENTER of ROTOR 

Op: CENTER of PRESSURE DISTRIBUTION 

1 : PHASE ANGLE of PRESSURE 
e .-ECCENTRICITY of SEAL 

APmax: MAXIMUM PRESSURE CHANGE 


Figure 10. - The assumption of the circumferential 
pressure distribution. 
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Figure 11. - An example of the analytic wave. 
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Figure 12. - Fluid force and phase angle of one chamber seal. 
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Figure 15. - Fluid force Figure 16. - Fluid force 
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Figure 23. - Damping coefficients in the 2nd chamber 
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Figure 24. - Spring coefficients in the 3rd chamber 
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Figure 25 


Damping coefficients in the 3rd chamber 
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Figure 26.- Spring coefficients in the 4th chamber 
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Figure 27. - Damping coefficients in the 4th chamber 
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